Complexity of Random Energy Landscapes, Glass Transition 
and Absolute Value of Spectral Determinant of Random Matrices 
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Finding the mean of the total number Ntot of stationary 
points for N— dimensional random energy landscapes is re- 
duced to averaging the absolute value of characteristic poly- 
nomial of the corresponding Hessian. For any finite N we 
provide the exact solution to the problem for a class of land- 
scapes corresponding to the "toy model" of manifolds in ran- 
dom environment. For N 1 our asymptotic analysis re- 
veals a phase transition at some critical value /i c of a control 
parameter fi from a phase with finite landscape complexity: 
Ntot ~ e NS , < jtt c ) > to the phase with vanishing 

complexity: E(/x > (x c ) = 0. Finally, we discuss a method of 
dealing with the modulus of the spectral determinant appli- 
cable to a broad class of problems. 



Characterising geometry of a complicated landscape 
described by a random function TL of N real variables 
x = (x±, xjv) is an important problem motivated by 
numerous applications in physics, image processing and 
other fields of applied mathematics [1]. The simplest, yet 
non-trivial task [2-6] is to find the mean number of all 
stationary points of 7i (minima, maxima and saddles) in 
a given domain of the Euclidean space by investigating 
the simultaneous stationarity conditions dkH = for all 
k = 1,...,N, with dk standing for the partial derivative 
-7T-. In this context the function 
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given by the sum of a purely deterministic quadratic piece 
characterized by a non-negative parameter /x > and 
of a random Gaussian function V(x) attracted consider- 
able interest for several independent reasons. For small 
N = 1,2 statistics of stationary points of (1) were in- 
vestigated long ago in a classical study of specular light 
reflection from a random sea surface [2] and addressed 
several times ever since in various physical contexts, see 
[3,4]. Most frequently one assumes Gaussian part to be 
isotropic with zero mean and correlations depending only 
on the Euclidean distance |xi — x 2 | 2 = Y^k=i( x ^k — X2k) 2 
and given by 



(V(xi)V(x 2 )) =Nf 
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with the brackets standing for the ensemble average. 



Recently much interest in Eq.(l) was boosted by re- 
interpreting it as the energy functional of a certain " toy" 
models describing elastic manifolds propagating in a ran- 
dom potential, see [7,8] and references therein. This type 
of models is known to display a very nontrivial glassy 
behaviour at low enough temperatures - an unusual off- 
equilibrium relaxation dynamics attributed to a complex 
structure of their energy landscape. Although particular 
dynamical as well as statical properties may differ sub- 
stantially for different functions f(x) (e.g. "long-ranged" 
vs. "short-range" correlated potentials, see [8]), the very 
fact of glassy relaxation is common to all of them. In fact, 
the same model admits an alternative interpretation as 
a spin-glass, with Xi being looked at as "soft spins" in a 
quadratic well interacting via a random potential V [7]. 
From this point of view it is most interesting to concen- 
trate on the limit of large number of "spins" N ^> 1. The 
experience accumulated from working with various types 
of spin-glass models [5] suggests that for the energy land- 
scape to be complex enough to induce a glassy behaviour 
the total number of stationary points Ntot (/-0 should grow 
exponentially with N as iVt t(/i) ~ exp7V£(/Lt). The 
quantity > in such a context is natural to call 

the landscape complexity. On the other hand, it is com- 
pletely clear that the number of stationary points should 
tend to N to t = 1 for very large /x when the random part 
is negligible in comparison with the deterministic one. 
In fact, when N — > oo we will find that a kind of sharp 
transition to the phase with vanishing complexity occurs 
at some finite critical value fi C) so that = as long 
as /i > /ic, whereas > for /i < fj, c and tends to 

zero quadratically when fi — > /i c . Such a transition is just 
the glass transition observed earlier in a framework of a 
different approach in [7,8]. 

We start with writing the number of stationary points 
of H in any spatial domain D as N^) = J D p(x) d^x, 
with jo(x) being the corresponding density of the sta- 
tionary points. The ensemble-averaged value of such a 
density can be found as 

P»„(x) =^|det (/z<5 fcl ik3 + dl M V) | Yl 6(fix k + d k V)^ 

where 5{x) and 5 mn stand for the Dirac's 5— function 
and Kroneker symbol, respectively. 

To evaluate the ensemble average we notice that 
for the Gaussian potential V the first derivatives d k V 
are Gaussian-distributed and are locally statistically 
independent of the second derivatives. Representing 
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the S— functions as Fourier integrals and exploiting 
{d m Vd n V) = a 2 5 mn , a 2 = — / (0) one can easily per- 
form the corresponding part of the averaging and arrive 
at: 
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|det(/iJ fcl , fc2 +H klk2 \) (3) 



where we introduced the matrix of second derivatives 
of the potential H klk2 = d\ x k V. Further changing 
H — > —H we see that the problem basically amounts to 
evaluating the ensemble average of the absolute value of 
the characteristic polynomial det(^/Ar — H) (a.k.a. spec- 
tral determinant) of a particular random matrix H . In 
particular, the total number of stationary points in the 
whole space is given by: 
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(\det{fiI N - H)\) 
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Whenever the physical problem necessitates to deal 
with the absolute value of the determinant, its presence 
considered to be a serious technical challenge, see [9] and 
references therein. In particular, an intensive work and 
controversy persists in calculating the so-called thermo- 
dynamic complexity of the free energy for the standard 
Sherrington-Kirkpatrick model of spin glasses [6] or its 
generalisations [5]. Several heuristic schemes based on 
various versions of the replica trick were proposed in the 
literature recently to deal with the problem, see discus- 
sion and further references in [6]. Despite some impor- 
tant insights, the present status of the methods is not yet 
completely satisfactory. 

In the present paper we propose two different methods 
of dealing with the modulus of determinants, both free 
from any mathematical uncertainty. The first method 
is specific for the problem in hand and is heavily based 
on the isotropy of the correlation function of the random 
field V in eq.(2). Exploitation of this fact provides one 
with a possibility to employ, after some manipulations, 
the standard methods of the random matrix theory [10] 
and find the explicit expression for the number N to t of 
stationary points for any spatial dimension N. 

To follow such a route we notice that the statisti- 
cal properties of the potential V result in the follow- 
ing second-order moments of the entries = 
1,...,N}: 



(HuHjm) = — [5ij5l m + S im Slj + SuSjm] (5) 

where we denoted J 2 = f (0). This allows one to write 
down the density of the joint probability distribution 
(JPD) of the matrix H explicitly as 
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where dH = Yii<i<j<N dHij and the proportionality 
constant can be easily found from the normalisation 
condition and will be specified later on. It is evident 
that such a JPD is invariant with respect to rotations 
H -» 0~ x HO by orthogonal matrices O G O(N), but it 
is nevertheless different from the standard one typical for 
the so-called Gaussian Orthogonal Ensemble (GOE) [10]. 
However, introducing one extra Gaussian integration it is 
in fact straightforward to relate averaging over the JPD 
(6) to that over the standard GOE. In particular, 



(|det {fiI N -H)\) = 

r°° dt 

J — o 



(7) 



e - w -(|dct[( A1 + Jt)I N -H }\) 



GOE 



where the averaging over Hq is performed with the GOE- 
typc measure: dHo Cm exp { — jj^TtHq }, with C/v = 
N 1 / 2 /[(2TTj 2 /N) N( - N + 1 ^ i 2 N ^ 2 } being the relevant nor- 
malisation constant. 

To evaluate the ensemble averaging in (7) in the most 
economic way one can exploit explicitly the mentioned 
rotational O(N)— invariance and at the first step in a 
standard way [10] reduce the ensemble averaging to the 
integration over eigenvalues Ai, X n of the matrix Hq. 
After a convenient rescaling Ai — ► J^2/N\i the resulting 
expression acquires the following form: 

/OO />OC 
d\x . . . / d\ N 
-oo J — OO 

(8) 
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x \{ \Xi - Xj I Y[ I ^/N/2(m + t) - Xi\e 
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where we denoted m — /i/J. One may notice that the 
above TV— fold integral can be further rewritten as a iV+1 
fold integral: 



e f (m+t) 



N+l 

N+i I I e 
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/ dXi... dX N+ i Yl 

J — oo J — 

x S (^Wj2(m + t) - A^+i) J] \ X i - \ 

i<j 

Such a representation makes it immediately evident that, 
in fact, the expectation value of the modulus of the de- 
terminant in question is simply proportional to the mean 
spectral density Vfq+i[m + t] (a.k.a one-point correlation 

function R[ N+l) y/N/2(m + t) , see [10]) of the same 

GOE matrix H Q but of enhanced size (N + 1) x (N + 1): 



{\det[(ti + Jt)I N -H ]\) GOE 
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(m+t)' 



VN+i[{m + 1)] 



vn[X] = — (TtS(XI n - H Q )) GOE 



(9) 



The last relation provides the complete solution of our 
original problem for any value of N, since the one-point 
function r[ N+1 ^ [x] is known in a closed form [10] for any 
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value of N in terms of the Hermite polynomials Hk (x) . 
In particular, for any odd integer N we have: 



|dct[(// + Jt)I N -H }\) 
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e 2 Hn + i (u)sign(a; — u) du 



(10) 



where we denoted x = y -y (m + t) for brevity For even 

integer iV one more term arises, see [10]. 

Being interested mainly in extracting the complexity 
= linijv^oo A^ 1 In N tot (p) we have to perform an 
asymptotic analysis of equations (4), (7), (9)-(10). In 
principle, one can employ the known large— N asymp- 
totics of the Hermite polynomials, but we find it more 
convenient to use an alternative, well-known representa- 
tion for the mean eigenvalue density: 



vn[X] = tt— Im-^ 



det (XI N - H ) 



Ntt dX b \ dct (X b I N - H ) 



U„=A 



(11) 



The ensemble average of the ratio of the two determi- 
nants can be easily found in the framework of the su- 
persymmetric approach [11]. Following a variant of this 
method we use for our analysis the following integral rep- 
resentation for the derivative of the ratio of two determi- 
nants featuring in Eq.(ll), see Eq.(46) in Ref. [12]: 

v N [X] oc Re ^ ^ e - NC ^G N (q, X) (12) 



where C(q) = \ + iXq — In (q) and 



G N (q) 
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dp 



(pi - q) 

3/2 
Pi 
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(P2 - q) , 



exp{--£( Pl )} (13) 
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d Pi — 3/2— l-Pi -P2|exp{-— £{p 2 )} 

P2 

The form of the above expressions Eqs. (12-13) suggests 
that the large- N asymptotics should be given by a saddle- 
point contribution in all integration variables. It turns 
out however that the situation is not that simple. To 
perform the asymptotic analysis accurately it is conve- 
nient first to get rid of the non-analyticity in the inte- 
grand by passing in the last expression to new variables 
< r < oo, — oo < 9 < oo by pi — re e ,p 2 = re~ e . 
Further introducing u = r (cosh 9 — 1) we arrive at 

f°° dr 

G N (q)<x / — cxp{-7V£(r)} (14) 
Jo r 

x / du(r + u)[(r-qf - 2qu}e' N ^+ 2u( - r+ ^ 
Jo 

The saddle-point value in u— variable is obviously u s — 
— (r + iX/2), but it can yield no contribution as long as 



Rer > along the contour of integration. The anal- 
ysis reveals that this is indeed the case for |A| < 2. 
Under that condition the u— integral is dominated by 
vicinity of its end point u = 0, whereas integrals over 
r and q are instead saddle-point domina ted, th e domi- 
nant saddle-points being r s = (— i\ + \/4 — A 2 )/2 and 
q s = (— i\ — V4 — A 2 )/2. Calculating the corresponding 
contribution we arrive, as expected, to the standard semi- 
circular spectral density i/[X] — ^\/4 — A 2 . If, however, 
the parameter A is such that |A| > 2, the situation turns 
out to be very differ ent. In that case both saddle-point 
values r s = i(— A± VA 2 — 4)/2 are purely imaginary, ne- 
cessitating a part of the steepest descent contour to be 
chosen along the imaginary axis Rer = 0. As a result, 
an additional contribution from the saddle-point u s turns 
out to be operative. Although such contribution is expo- 
nentially small in comparison with one dominated by the 
vicinity of u = 0, it is the only one which survives after 
taking Re in (12). Taking into account the saddle-pont u s 
induces modifications of the relevant exponential term in 
r— variable, which now becomes expjA^^ + lnr — ^-]} 
and replaces the former expression exp {— N£(r)}. The 
relevant saddle-point for r then turns out to be r s = —i 
as long as A > 2 and it results in exponentially small 
("instanton") value for the spectral density: 



v[X] oc exp ■ 
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tAVA 



2 - 4 - In 



X + VX 2 



A > 2 



(15) 



where we only kept factors relevant for calculating the 
complexity in the limit of large N . 

Finally, we employ the relation (9) between the mean 
spectral density and the expectation value of the modulus 
of the spectral determinant for GOE matrices, and sub- 
stitute the resulting expression into the integral (7). In 
the latter we can again exploit the saddle-point method 
for asymptotic analysis. For < m < 1 the relevant 
saddle point is t s — m satisfying 0<A s =t s +m<2, 
and validati ng the u se of the semicircular spectral density 
v[X s ] = ^ y/A — A 2 in the calculation. This yields 

(| det (fj,I N -H)\) oce^^Vl-rrt 2 , 0<m<l (16) 

For m > 1, however, it turns out that one has to use 
Eq.(15) for the spectral density. The corresponding 
saddle-point value t s in the t— integral is given by the 
solution of the equation m = ^(X s + \J A 2 — 4) for the 
variable X s = t s + m. The solution is easily found to be 
simply X s = m + mT 1 (note that A s > 2 ensuring consis- 
tency of the procedure) which yields the resulting value 
for the modulus of the determinant to be given by 



det (/i/jv — H)\) oc e 



Nlnm 



m > 1 



(17) 



Invoking our basic relation Eq.(4) for N to t we see that 
the landscape complexity S(/i) of the random potential 
function (1) is given by: 
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X(ji) = l (j2-l) - T). ,, „, = ./ (18) 



(19) 



Earlier works referred to the critical value fi c — J as, 
on one hand, signalling the onset of a nontrivial glassy 
dynamics [8], and, on the other hand, corresponding to 
the point of a breakdown of the replica-symmetric solu- 
tion [7]. Our calculation provides an independent sup- 
port to the point of view attributing both phenomena 
to extensive number of stationary points in the energy 
landscape. At the critical value the complexity vanishes 
quadratically: — > /j, c ) oc (yu c — fi) 2 /fi%- 

Finally, let us very shortly discuss an alternative, less 
model-specific technique of evaluating the absolute value 
of the spectral determinant. It is based on the following 
useful identity , see e.g. [13]: 



| det (fj,I N -H)\ = 



(20) 



lim ■ 



[det {^I N -H)f 



y/dct (Jp - ie)I N - H)y/det (Jp + ie)I N - H) 

valid for any matrix H with purely real eigenvalues. For 
the particular case of real symmetric matrices H one can 
represent the two factors in the denominator of the right 
hand side in terms of the gaussian integrals absolutely 
convergent as long as e > 0. Further representing the 
determinantal factors in the numerator in terms of the 
Gaussian integral over anticommuting (Grassmann) vari- 
ables we thus get a bona fide supersymmetric [11] object 
to be analysed. Simultaneous presence in the starting 
expression both /i* = fi ± ie and fj, makes the calcula- 
tion in this case more involved in comparison with just 
a simple ratio of two determinats, as in (11). It is nev- 
ertheless an important fact that possibility to perform 
the ensemble average explicitly exists whenever matrix 
entries of H are Gaussian-distributed, not requiring any 
matrix invariance or even independence of the matrix 
entries. Similar strategy may be even employed when H 
is a stochastic differential operator with certain Gaussian 
part, as in the notorously difficult case of the random field 
Ising model [14]. For this reason it is natural to hope that 
the suggested method may appear of certain utility be- 
yond the present model, e.g. when discussing free energy 
landscapes for spin-glass related problems [5,6]. 

Full account of the procedure will be presented in a 
separate publication and here we just quote the result 
of the ensemble averaging for any finite value of N (cf. 
Eq.(12)): 



| det {pI N -H )\) GOE oc 



(21) 



>< r r P^^-^f^ N[c(qi)+c[q2)] 

J-ooJ-oo (<H<12) 2 

x / dpi(pi + qi){pi + q 2 ) / dp 2 (p 2 - qi)(P2 - 92) 
Jo Jo 

N 



x [^3/2 K [Ne( Pl +P2)]cx P -y a Pl ) + a P2 ) 



where m = fi/J, Kq(z) stands for the Macdonald func- 

2 

tion, C(q) = 3j- + imq — In (q), and the bar stands for 
the complex conjugation. Note, that because of the log- 
arithmic divergency of the Macdonald function at small 
arguments one should keep n — Ne finite when perform- 
ing an accurate asymptotic analysis of this integral, and 
only then set k — 0. The resulting expression reproduces 
all the features found by us earlier in the present paper. 

In summary, we calculated the mean total number iV tot 
of stationary points for a N— dimensional potential con- 
sisting of a quadratic well of strength \i and of a ran- 
dom gaussian piece V. In particular, for ./V — > 00 we 
found that the potential is characterised by finite land- 
scape complexity: N to t ~ e NTl , S > as long as [i < fi c , 
and for /j, — ► ji c the complexity X vanishes quadratically. 
Finally, we discuss a general method of calculating the 
mean absolute value of the spectral determinant. 
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